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ABSTRACT. 

lAU 2000 resolutions on the reference frames set up a solid theoretical foundation for implementing 
general relativity in astronomical data processing algorithms and for unambiguous interpretation of mea- 
sured relativistic effects. We discuss possible directions for further theoretical development of the lAU 
resolutions aimed to take into account the decadal progress in observational techniques and computer- 
based technologies. We address the following subjects: 1) space-time transformations and the structure 
of the metric tensor; -2) PPN parameters and gauge invariance of equations of motion; -3) astronomical 
reference frames for cosmological applications. 

1 INTRODUCTION 

Experimental exploration of the nature of space-time demands establishment of a common theoretical 
platform linking a theory of gravitational field to astronomical observations. This platform should be 
build on the basis of a complete theory of gravity like general theory of relativity, that describes both 
the properties of space-time, gravitational field and observables. New generation of microarcsecond 
astrometry satellites like SIM Lite ^ and a cornerstone mission of ESA - Gaia, require such a novel 
approach for an unambiguous interpretation of astrometric data obtained from the on-board optical 
instruments. Advanced inertial reference frame is required for unambiguous physical interpretation of 
gravitomagnetic precession of LAGEOS satellite and LLR observations [1]. Recent breakthroughs in 
technology of drag-free satellites, clocks, lasers, optical and radio interferometers and new demands 
of experimental gravitational physics [2, 3] make it necessary to incorporate the parameterized post- 
Newtonian formalism [4] to the procedure of construction of relativistic local frames around Earth and 
other bodies of the solar system [5, 6]. The domain of applicability of the lAU relativistic theory of 
reference frames [7] is to be also extended outside the solar system [8] to take into account the impact of 
the Hubble expansion on the solutions of the gravity field equations and the equations of motion of the 
bodies. 

In what follows, Latin indices takes values 1,2,3; the Greek indices run from to 3. Repeated indices 
imply the Einstein summation rule. The unit matrix Sij = diag(l,l,l) and the fully anti-symmetric 
symbol eijk is subject to £123 = 1. The Minkowski metric ijap = diag(— 1, 1, 1, 1). Greek indices are raised 
and lowered with the Minkowski metric, Latin indices are raised and lowered with the unit matrix. Bold 
italic letters a, b, etc., denote spatial vectors. A dot and a cross between two spatial vectors denote the 
Euclidean scalar and vector products respectively. Partial derivative with respect to spatial coordinates 
are denoted as d/dx^ or V. 

2 STANDARD lAU FRAMEWORK 

The lAU resolutions arc based on the first post-Newtonian approximation of general relativity which is 
a conceptual basis of the fundamental astronomy in the solar system [9] . Barycentric Celestial Reference 



^Thc Astro2010 Decadal Survey (available at http://sites.nationalacadeinies.org/bpa/BPA_049810) did not recom- 
mend SIM Lite for development this decade. 



System (BCRS), = (ct, x), is defined in terms of a metric tensor ga/s with components 



500 = + — - — + O c-5 , 1 

50. = -^ + 0(c-5), (2) 

= % (^1 + ^) + 0(c-4) . (3) 

Here, the post-Newtonian gravitational potentials w and are defined by solving the gravity field 
equations 

Dw = -inGa , (4) 

Dw' = -AiiGa' , (5) 

where □ = -c~'^d'^/dt^ + is the wave operator, a = c^^j-j^oo _^ yss^^ ^, ^ c-^T"*, and T^''' are the 
components of the stress-energy tensor of the solar system bodies, T^'^ — T^^ + T^^ -t- T^^ . 
Equations (4), (5) are solved by iterations 

.'(t,.) ^ G/4^i^ + 0(c-), (7) 



which arc to be substituted to the metric tensor (l)-(3). Each of the potentials, w and w' , can be linearly 
decomposed in two pieces 

w = We + w , (8) 
= w% + w\ (9) 

where we and w'^ are BCRS potentials depending on the distribution of mass and current only inside 
the Earth, and we and w^^ are gravitational potentials of external bodies. 

Geocentric Celestial Reference System (GCRS) is denoted X"' = {cT,X). It has the metric tensor 
Gai3 with components 

Goo = -l + ^--^ + 0(c-5) , (10) 

Go, = -^ + 0(c-5), (11) 
( 2W\ 

G,, = (^1 + — j + 0(c-4) . (12) 

Here W = VF(r, X) is the post-Newtonian gravitational potential and VF*(r, X) is a vector-potential 
both expressed in the geocentric coordinates. They satisfy to the same type of the wave equations (4), 
(5). Planetocentric metric for any planet can be introduced in the same way as the GCRS. 
The geocentric potentials. We and W^|;, arc split into three parts 

w{T,x) = VK£(r,x) + VKkin(r,x) + Wdy„(r,x) , (13) 

W\T,X) = Wk{T,X) + WUT,X) + W^^,,{T,X) . (14) 

associated respectively with the gravitational field of the Earth, external tidal field and kinematic inertial 
force. lAU resolutions imply that the external and kinematic parts must vanish at the geocenter and 
admit an expansion in powers of X [7]. Geopotentials We and W^ are defined in the same way as we and 

in equations (6)-(7) but with a and cr* calculated in the GCRS. They are related to the barycentric 
gravitational potentials we and w''^ by the post-Newtonian transformations [10, 7]. 

The kinematic contributions are linear in the GCRS spatial coordinates X 

M^kin = Q^X^ , Wi,^ = i ch.,,{nP - QP,,,) X^ , (15) 



where Qi characterizes a deviation of the actual world line of the geocenter from a fiducial world line of 
a hypothetical spherically-symmetric Earth [11] 

Q, = OM^e) - + 0(c-2) . (16) 

Here = dv'^^ /dt is the barycentric acceleration of the geocenter. Function rip^gc describes the relativistic 
precession of dynamically non-rotating spatial axes of GCRS with respect to reference quasars 

f^prcc = ^£.jfe [-lv'Edkw{XE) + 2dkW^XE)-^v'EQ''y (17) 

The three terms on the right-hand side of this equation represent the geodetic, Lense-Thirring, and 
Thomas precessions, respectively [11, 7]. Dynamic potentials Wdyn and W^^^ are generalizations of the 
Newtonian tidal potential in the form of a polynomial starting from the quadratic with respect to X 
terms. 



3 lAU SCALING RULES AND THE METRIC TENSOR 

The coordinate transformations between the BCRS and GCRS arc found by matching the BCRS and 
GCRS metric tensors in the vicinity of the world line of the Earth by making use of their tensor properties. 
The transformations are written as [11, 7] 

T = t-\[A + VE-rE] + \\B + BW'E + B''r'Er%], (18) 



X' = + 1 



c 

where te — x — xe-, functions A, B, B*, i?*^ obey equations 
dA 1 



1:VeVe ■ Ve + w{xE)rE + rEttE ■ rE - \a:Er% 



(19) 



^ = ^^E+^ixE). (20) 
^ = -\v%~^vIw{xe)+Av'ew' + ^w^{xe), (21) 
B' = -~]^iiiv\, + Aw'{xE)-iv'Ew{xE), (22) 



B'^ = -v'EQj+2djw\xE)-v'Ed,wixE) + -S'^7jj{xE), (23) 



where x^, v^^, and are the BCRS position, velocity and acceleration vectors of the Earth, the ovcrdot 
stands for the total derivative with respect to t, and one has neglected all terms of the order 0{r^). 
Earth's orbit in BCRS is almost circular. This makes the right side of equation (20) looks like 

-Ve + w{xe) = c^Lc + (periodic terms) , (24) 

where the constant Lc and the periodic terms have been calculated with a great precision in [13]. For 
practical reason the lAU 2000 resolutions recommend to re-scale the BCRS time coordinate t to remove 
the constant Lc from equation (24). The new time scale is called TDB, and it is defined by equation 

tTDB^t{l-LB) , (25) 

where a constant Lb = Lc + ^c is used, instead of Lc, in order to take into account the additional linear 
drift Ac between the GCRS time T and the proper time of clocks on gcoid, as explained in [12, 13]. Time 
re-scaling changes the Newtonian equations of motion of planets and light. In order to keep the equations 
of motion invariant entails re-scaling of spatial coordinates and masses of the solar system bodies. These 
scaling transformations are included to lAU 2000 resolutions [7]. However, the re-scaling of masses, times 
and spatial coordinates affects the units of their measurement - the procedure that led to a controversial 
discussion [14]. The change in units can be avoided if one looks at the scaling laws from the point of view 
of transformation of the metric tensor. 



The thing is that the scahng of time and space coordinates can be viewed as a particular choice of 
the GCRS metric tensor. Indeed, equation (15) is a solution of the Laplace equation which is defined up 
to an arbitrary function of time Q = Q{t) that can be incorporated to 



dA{t) 
dt 



= Q 



1 2 
2^^ 



w{xe) 



(26) 
(27) 



and, if one chooses Q = c^Lc, it eliminates the secular drift between times T and t without explicit 
re-scaling of the time t, which is always measured in SI units. It turns out that Blanchet-Damour [15] 
relativistic definition of mass depends on function Q and is re-scaled automatically in such a way that the 
Newtonian equations of motion remain invariant [6] . Introduction of Q to function Wkin appropriately 
transforms the gij component of the GCRS metric tensor that is formally equivalent to the previously- 
used re-scaling of the GCRS spatial coordinates. One concludes that introducing the function Q to the 
GSRC metric tensor without apparent rc-scaling of coordinates and masses can be more preferable in 
updated version of the lAU resolutions as it allows us to keep the SI system of units without changing 
coordinates and masses. Similar procedure can be developed for the topoccntric metric tensor to take 
into account the linear drift existing between GCRS time T and the atomic clocks on geoid [12, 13]. 

4 PARAMETERIZED COORDINATE TRANSFORMATIONS 



The parameterized post-Newtonian (PPN) formalism [4] is not consistent with the lAU resolutions. It 
limits applicability of the resolutions in testing gravity theories. PPN equations of motion depend on two 
parameters, /3 and 7 [16] and they are presently compatible with the lAU resolutions only in the case of 
/3 = 7 = 1. Rapidly growing precision of astronomical observations as well as advent of gravitational- 
wave detectors urgently demand a PPN theory of relativistic transformations between the local and global 
coordinate systems. 

PPN parameters /3 and 7 arc characteristics of a scalar field which makes the metric tensor different 
from general relativity. In order to extend the lAU 2000 resolutions to PPN formalism one used a general 
class of Brans-Dicke theories [17] based on the metric tensor gap and a scalar field (j) that couples with 
the metric tensor via function 0{(f>). Both and 0{(f>) arc analytic functions which can be expanded in a 
Taylor scries about their background values and 9. 

The parameterized theory of relativistic reference frames in the solar system is built in accordance to 
the same rules as used in the lAU resolutions. The entire procedure is described in papers [5, 6]. The 
PPN transformations between BCRS and GCRS are found by matching the BCRS and GCRS metric 
tensors and the scalar field in the vicinity of the world line of the Earth. They have the following form 
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where te = x — xe, and functions A{t), B{t), B^{t), B'^^{t) obey 
1 



dA 

m 

dB 

B' 



1 , 



- Q{xe), 

7 + i j w| w{xe) + 2(1 H- 7) „^ + f /3 - 1 ) w^{xe), 



vl v'e + 2(1 + 7) w'{xe) - (1 + 27) v'e w{xe), 



-VeQj + (1 + -/)djw''{xE) - -fVEdjW{xE) + - S'^wixE) 



(28) 
(29) 

(30) 
(31) 
(32) 
(33) 



These transformations depends explicitly on the PPN parameters /3 and 7 and the scaling function Q, 
and should be compared with those (18)-(23) currently adopted in the lAU resolutions. 

PPN parameters /3 and 7 have a fundamental physical meaning in the scalar-tensor theory of gravity 
along with the universal gravitational constant G and the fundamental speed c. It means that if the 
parameterized transformations (28)- (33) are adopted by the lAU, the parameters /3 and 7 are to be 
considered as new astronomical constants which values have to be determined experimentally. 



5 MATCHING lAU RESOLUTIONS WITH COSMOLOGY 



BCRS assumes that the solar system is isolated and space-time is asymptotically flat. This idealization 
will not work at some level of accuracy of astronomical observations because the cosmological metric 
has non-zero Riemannian curvature [18]. It may turn out that some, yet unexplained anomalies in the 
orbital motion of the solar system bodies are indeed associated with the cosmological expansion [19]. 
Moreover, astronomical observations of cosmic microwave background radiation and other cosmological 
effects requires clear understanding of how the solar system is embedded to the cosmological model. 
Therefore, it seems reasonable to incorporate the cosmological metric to the lAU resolutions. 

The gravitational field of the solar system has to approximate the cosmological metric tensor at 
infinity, not a flat metric. The cosmological metric has a number of parameters depending on visible 
and dark matter and on the dark energy. One considered a universe, driven by a scalar field imitating 
the dark energy 0, and having a spatial curvature equal to zero [20, 21]. The universe is perturbed by a 
localized distribution of matter of the solar system. The perturbed metric tensor reads 

gap = a^ir])fal3 , fap = Va/S + , (34) 

where the perturbation hais of the background metric gap = a^rjap is caused by matter of the solar 
system, 0(77) is a 'radius' of the universe depending on the conformal time 77 related to coordinate time t 
by simple differential equation dt — a(ri)dri. A linear combination of the metric perturbations 



7 

where h = rj^^hap, is more convenient for calculations. 
One imposes a cosmological gauge given by [20, 21] 

r^ip ^ 2H^SS , (36) 

where a vertical bar denotes a covariant derivative with respect to the background metric gap, ^ = </>/a^, 
H = a/a is the Hubble parameter, and the overdot denotes a time derivative with respect to time rj. 
The gauge (36) generalizes the harmonic gauge of the lAU resolutions for the case of the expanding 
universe. The gauge (36) drastically simplifies the field equations. Introducing notations 700 = 4z/;/c^, 
7oi = —Aw'^/c^, and jij = 4i/;*-'/c^, and splitting Einstein's equations in components, yield 

nx-2Hd^X+lH\ = -AirGa, (37) 

Dw - 2Hd^w = -AnGa - iH^x , (38) 

Dw' - 2Hd^w' + H^w' ^ -AnGa' , (39) 
Dw'^ - 2Hdr,w'^ = -AttGT'^ , (40) 

where drj = d/drj, □ = — c~^9^ -I- V^, x = — 1^/2, the Hubble parameter H = a/a = 2/77, densities 
a = c~^(T°° -I- T*'*), cr* = c~^T^^ with T"^ being the tensor of energy- momentum of matter of the solar 
system defined with respect to the metric fap- These equations extend the equations (4), (5) of the lAU 
resolutions to the case of expanding universe. 

Equation (37) describes evolution of the scalar field while equation (38) describes evolution of the 
scalar perturbation w of the metric tensor. Equation (39) yields evolution of vector perturbations of the 
metric tensor, and equation (40) describes TT gravitational waves emitted by the solar system. Equations 
(37)-(40) depend on the Hubble parameter and can be solved analytically. The Green functions for these 
equations have been found in [20, 21] and solutions can be smoothly matched with the BCRS metric (6), 
(7) of the lAU resolutions. 

This work was promoted by the Research Council Grant FIT-11-020 of the University of Missouri- 
Columbia. The author is thankful to N. Capitaine and the LOC of Journees 2010 for travel support. 
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